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\^ Abstract 

O; 

^J^ I We formulate the hadronic tensor VF^jy of deep inelastic scattering in the 

,__! ' path-integral formalism. It is shown that there are 3 gauge invariant and topo- 

Q\ ■ logically distinct contributions. Besides the valence contribution, there are two 

^^ , sources for the sea ~ one in the connected insertion and the other in the dis- 

r^ I connected insertion. The operator product expansion is carried out in this 

II ■ formalism. The operator rescaling and mixing reveal that the connected sea 

Oh. partons evolve the same way as the valence, i.e. their evolution is decoupled 

(— I I from the disconnected sea and the gluon distribution functions. We explore the 

L| ' phenomenological consequences of this classification in terms of the small x be- 

•rH , havior, Gottfried sum rule violation, and the flavor dependence. In particular, 

rS I we point out that in the nucleon u and d partons have both the connected and 

C^ ' disconnected sea contributions; whereas, s parton has only the disconnected 

sea contribution. This difference between u + d and s, as far as we know, has 

not been taken into account in the fitting of parton distribution functions to 

experiments. 
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1 Introduction 

In the last decade, the surprising results of a small quark spin content (flavor-singlet 
9a) mi ^^^ ^^^ discovery that u ^ d in the nucleon from deep inelastic scattering 
have focused people's attention on the interplay between the parton model at high 
energies and the hadronic structure at low energies. The connection is often made 
through the operator product expansion which relates the sum rules of parton dis- 
tribution functions to the forward matrix elements. The latter can be obtained from 
low energy experiments. 

In the parton model, the dynamical quark degrees of freedom are taken to be the 
valence and the sea. Whereas, in the hadronic models, the degrees of freedom involve 
the valence and the meson cloud. The classification of these dynamical degrees of 
freedom in deep inelastic scattering has been made in the path integral formalism 0, 
Q . It is revealed that there are two sources for the sea quark. One is in the connected 
insertion and the other in the disconnected insertion. Their relations to the meson 
cloud in the hadronic models for hadrons near the rest frame and chiral symmetry 
have been clarified and extensively explored in the context of hadronic models in terms 
of the form factors, hadron masses, and matrix elements, i.e. low-energy quantities 
which can be calculated in the two- and three-point functions. In addition, it is 
shown that when both the connected sea (referred to as cloud in Ref.jQ) and the 
disconnected sea quarks are eliminated in a valence QCD theory, the valence quark 
picture with SU(6) symmetry emerges. In this paper, we shall derive the operator 
product expansion in the path-integral formalism and explore the phenomenological 
consequences of this classification of the parton degrees of freedom. Section 2 is 
on the path-integral formalism of the hadronic tensor and the classification of the 
parton degrees of freedom. Section 3 shows how to carry out the operator product 
expansion in the path-integral formalism and to derive parton evolution equations 
through operator rescaling and mixing. Section 4 explores the phenomenological 
consequences. We shall show that the small x behavior of the connected sea is different 
from that of the disconnected sea and the violation of the Gottfried sum rule, i.e. 
u ^ d, comes only from the connected sea at the flavor SU{2) limit. Finally, we 
emphasis that u+d in the nucleon has both the connected sea and the disconnected sea 
contributions; whereas s has only the disconnected sea contribution. This difference 
has not been parametrized in extracting the parton distribution functions from the 
experiments. The conclusion is given in Section 5. 



2 Path-Integral Formalism 

The deep inelastic scattering of a muon on a nucleon involves the hadronic tensor 
which, being an inclusive reaction, involves all intermediate states 
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W,M\^) = 7TJ^Y.i'^^f^\Pn-p-q){N\J,mn){n\JM\N) 
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Since deep inelastic scattering measures the absorptive part of the Compton scatter- 
ing, it is the imaginary part of the forward amplitude and can be expressed as the 
current- current correlation function in the nucleon, i.e. 






It has been shown p, |^ that the hadronic tensor W^y{q^, u) can be obtained from 
the Euclidean path-integral formalism where the various parton dynamical degrees 
of freedom are readily and explicitly revealed. In this case, one considers the ratio 
of the four-point function ^{ONit) J ff e-^«^-^J^(x, ^2)^^(0' ^i)C7v(0)) and the two- 
point function {0^(1 — (^2 — ii))O7v(0)), where 0^(1;) is an interpolation field for the 
nucleon with momentum p at Euclidean time t. 

As both t — t2 » 1/AEp and ti » 1/AEp, where AEp is the energy gap between 
the nucleon energy Ep and the next excitation (i. e. the threshold of a nucleon and 
a pion in the p-wave), the intermediate state contributions will be dominated by the 
nucleon with the Euclidean propagator e~^p(*~*^*2-ti)) Hence, 
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^^^ < N\ J ^e-^'^-^Ux,h)MO,h)\N >, (3) 

where t = t2 — ti and / is the transition matrix element (0|O7v|iV), and V is the 
3-volume. Inserting intermediate states, Wp_^{q'^,T) becomes 

ZMn n 

(4) 
To go back to the delta function S{En — Ep + u) in Eq. (P, one needs to carry out 
the inverse Laplace transform |^, ^ 

1 rc+ioo , 

W^,{q\ u) = - dre^^W^^iq', r), (5) 
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with c > 0. This is basically doing the anti-Wick rotation back to the Minkowski 
space. 

In the Euclidean path- integral formulation of Wf^^iQ"^, t) in Eq. (^), contributions 
to the four-point function can be classified according to different topologies of the 
quark paths between the source and the sink of the proton. They represent different 
ways the fields in the currents J^ and J^, contract with those in the nucleon interpo- 
lation operator On- Fig. 1(a) and 1(b) represent connected insertions (C.I.) of the 
currents. Here the quark fields from the interpolators On contract with the currents 
such that the quark lines flow continuously from t = to t = t. Fig. 1(c), on the other 
hand, represents a disconnected insertion (D.I.) where the quark fields from J^ and 
Ji/ self-contract and are hence disconnected from the quark paths between the pro- 
ton source and sink. Here, "disconnected" refers only to the quark lines. Of course, 
quarks dive in the background of the gauge field and all quark paths are ultimately 
connected through the gluon field. 
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Figure 1: Quark skeleton diagrams in the Euclidean path integral formalism for 
evaluating W^^ from the four-point function defined in Eq. (H). These include the 
lowest twist contributions to W^^j. (a) and (b) are connected insertions and (c) is a 
disconnected insertion. 



Fig. 1 represents the contributions of the class of "handbag" diagrams where 
the two currents are hooked on the same quark line. These contain leading twist 
contributions in deep inelastic scattering. Other contractions with the two currents 
hooking on different quark lines involve only higher twist operators and thus will 
be suppressed in the Bjorken limit [^]. They are shown in Fig. 2. We will neglect 
these "cat's ears" diagrams from now on. We should stress that these diagrams in 



Figs. 1 and 2 are not Feynman diagrams to repesent the forward Compton scattering 
amplitude T^y{q^, v) and should not be read as such. Rather, they are path- integral 
diagrams needed to formualte W^u{(f,T^) which is the imaginary part of T^yi^q^.v) 



(see Eq. 



or its s-channel discontinuity. 
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Figure 2: Quark skeleton diagrams similar to those in Fig. 1, except that the two 
current insertions are on different quark lines. They give higher twist contributions 
ioW, 
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In the deep inelastic limit, the Bjorken scaling implies that the current product 
(or commutator) is dominated by the light-cone singularity of a free-field theory, i. e. 
1/x^ where x^ ~ O^l/Q"^). Among the time- fixed diagrams in Fig.l, Fig. 1(a) /1(b) 
involves only a quark/antiquark propagator between the currents; whereas. Fig. 1(c) 
has both quark and antiquark propagators. Hence, there are two distinct classes of 
diagrams where the sea quarks contribute. One comes from the D.I. ; the other comes 
from the C.I.. It is usually assumed that connected insertions involve only "valence" 
quarks which are responsible for the baryon number. This is obviously not true, there 
are also quark- antiquark pairs in the C.I. To define the quark distribution functions 
more precisely, we shall call the quark/antiquark distribution from the D.I. (which 
are connected to the "valence" quark propagators and other quark loops through 
gluons) the "disconnected sea" quark/antiquark. We shall refer to the antiquark in 
the backward time going quark propagator between ti and ^2 in Fig- 1(b) as the 
"connected sea" antiquark. On the other hand, the quark in the time forward prop- 
agator between ^2 and ti in Fig. 1(a) includes both the valence and the "connected 
sea" quarks. This is because a quark propagator from t = to t = t(t > 0) involves 
both the time forward and backward zigzag motions so that one cannot tell if the 



quark propagator between ^2 and ti is due to the valence or the connected sea. All 
one knows is that it is a quark propagator. In other words, one needs to consider 
connected sea quarks in addition to the valence in order to account for the produc- 
tion of quark-antiquark pairs in a connected fashion (Fig. 1(a)); whereas, the pair 
production in a disconnected fashion is in Fig. 1(c). 

We should stress that this separation into three topologically distinct classes of 
path- integral diagrams is gauge invariant. Notice that all the quark propagators are 
sewed together in a trace over color. In a perturbative illustration of the distinction 
between Fig. 1(b) and Fig. 1(c), one may consider the time-ordered perturbation 
where Fig. 1(c) represents the vacuum polarization contribution as a disconnected 
insertion in a direct diagram. The corresponding exchange diagram where the quark 
in the loop in Fig. 1 (c) is exchanged with one in the "valence" will lead to a connected 
insertion which falls in the class of Fig. 1(b) |^, |^. However, the separation depends 
on the momentum frame of the nucleon, although the sum which corresponds to the 
full physical W^y{q^, v) does not. For example, when the quark/ ant iquark propagator 
between the currents is either from the nucleon interpolation field or pair-produced 
before the current J^ at ti, i.e. it is pre-existing in the wavefunction, then it is not 
suppressed in the large momentum frame. Whereas, if it is pair-produced by the 
current Jp, then it is suppressed by \p\^ where |p| is the momentum of the nucleon. 
This has been known since current algebra sum rules were studied at large |p| 0. 

Since the parton model acquires its natural interpolation in the large momen- 
tum frame of the nucleon, i.e. |p| > |g|, the parton distribution is then defined via 
i^W2{Q'^,i') — > F2{x,Q'^) = xY.i^fi.flii.x^Q'^) + Qi^XyQ"^)) in the large momentum 
frame. Here x is the Bjorken scaling variable x = Q'^ /2mv. Given the specific time- 
ordering in Fig.l, Fig. 1(a) /1(b) involves only a quark/antiquark propagator between 
the currents; whereas. Fig. 1(c) has both quark and antiquark propagators. Conse- 
quently, the parton densities for the u and d antiquarks in the nucleon come from two 
sources, i.e. for the case of u, 

m(x, Q^) = Ucs{x, Q^) + Udsix, Q^), (6) 

where Ucs{x,Q'^) is the "connected sea" (CS) u anti-parton distribution from the 
C. I. in Fig. 1(b) and Uds{x,Q'^) denotes the "disconnected sea" (DS) u anti-parton 
distribution from the D. I. in Fig. 1(c). Similarly, d has two components. The strange 
and charm partons, on the other hand, only appear in the D. I. in Fig. 1(c). Thus, 

six,Q') = Sdsix,Q^). (7) 

One can prove Eq. (^ this way. First, it is imperative to note that the haronic tensor 
W^^i^q"^, t) in Eq. (|^) does not depend on the specific form of the interpolation field, 
except its quantum numbers. In fact, the interpolation- field dependent transition 
matrix element / = {0\On\N) drops out in the ratio of the four-point to two-point 
functions in Eq. (^. As such, one can use the simplest interpretation field of the 



nucleoli which involves only the valence quark field. For example, ON{t) can be taken 
to be the two u and one d quark fields with nucleon quantum numbers, 



for the proton. Since the interpolation field On does not involve strange quarks, the 
strange parton contribution can only come from the vacuum polarization due to the 
external currents J^ and Jy, in other words, the D. I. in Fig. 1(c). As a corollary, 
one can also prove it for the case if one uses the O^ss as the interpolation field, 
for example. In this case, there are two classes of path-integral diagrams. One class 
involves the C 1. where the strange quark fields in the currents S7^s and s'-f^s contract 
with those strange quark fields in the interpolation field Onss for the nucleon source 
and sink. This class of diagrams does not project to the nucleon as its lowest mass 
state, since the physical states it projects to will involve 5 valence quarks, i.e. uudss. 
Instead, it will project to states such as the nucleon and a scalar ss meson. Since 
they all have masses higher than the nucleon, they will be exponentially suppressed 
relative to the nuelcon as the time separation t — t2 and ti in Eq. (^ are large. The 
other class involves a D. I. where the strange quark fields in the nucleon source and 
sink self contract, so are the strange quark fields in the currents. This will project to 
the nucleon state with uud as the valence quarks. Since the transition matrix element 
/ = {Q\Onss\N) is divided out in the ratio in Eq. (^, it yields the same result as 
that obtained with O^ as the interpolation field. Thus, in either case, the strange 
parton contribution comes only from the D. I. in Fig. 1 (c). 

Similarly, the u and d partons have 2 sources, i.e. 

m(x, Q"^) = u^,+cs{x, Q^) + Uds{x, Q^), (9) 

where Ui,+cs{x,Q^) denotes the valence and CS u partons and Uds{x,Q'^) denotes the 
DS u parton and they are from Fig. 1(a) and Fig. 1(c) respectively. Again these 
two components applies to the d partons; whereas, the s parton has only the DS 
component. 



3 Operator product expansion (OPE) 

In the Minkowski space, the operator product expansion (OPE) is carried out in 
the unphysical region of T^^, which is defined with the time-ordered product of the 
currents. How does one carry this out in the Euclidean path-integral formulation? It 
turns out that because W^^i^q"^ , r) is defined in the Euclidean path-integral (Eq. (|D), 
it requires several steps to get to T^^, in the Minkowski space. On the other hand, it 
is relatively easy to do so because it entails a simple Taylor expansion of functions as 
opposed to dealing with operators in the usual OPE, as we shall see. 
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Considering Fig. 1(a) first, the three-point function in Eq. (^ involves the fol- 
lowing expression 

W^,,{q\T) oc I d[A]detM{A)e-^^Tr[- ■ ■ M-\tM) j d^xe~"^-^ i-^^M-\t2,t2-T)i^,M-\t2-T,Q) ■ ■ ■], 

(10) 
where Sg is the action for the gluon field A, M ^ is the quark propagator with 
arguments labeled by the Euclidean time. The spatial indices are implicit and have 
been integrated over to give the nucleon a definite momentum \f)\ and a momentum 
transfer q. x and r are the spatial and time separations of the two currents J^ and Jy. 
The trace is over the color and spin indices. The expression in Eq. (|T^) exhibits the 
part of the result from the quark line on which the currents are attached. The other 
two quark propagators and the nucleon interpolation field operators are indicated by 
the dots. 

Similar to the usual OPE derivation |^, we shall consider the most singular part 
of the quark propagator between the currents in Fig. 1. In the DIS limit where 
both the momentum transfer |g| and energy transfer v —>■ oo, the leading singularity 
comes from the short- distance part in PV^iy(g^, r) where \x\ and r — »• 0. Therefore, we 
replace the quark propagator between the currents with the free massless propagator 

M-fe,«.-.)^jL_^. („, 

We also carry out the Taylor expansion of the propagator M~^{t2 — r, 0) for small r 
and X 

M-\t2 - r, 0) = e^^-^+"^-M-i(t2, 0), (12) 

where D is the covariant derivative. With these substitutions, the corresponding 
hadronic tensor Wfj,y{q^,v) from Fig. 1(a) after the Fourier transform in space and 
Laplace transform in r (Eq. (^)) is given as 

W^,{q\ v) oc Tr[ ■ ■ M-\t, ts) i^^ ^'''}^ ^ '^h {u+D,-\q+iD\) i-f,M-\t2, 0) ■ ■ ■]. 

\q + iD\ 

(13) 

Since W^y{q^,u) is the imaginary part of T^j,, i. e. W^y{q^,u) = -/mT^j,(g^, z/), 
one can use the dispersion relation to obtain T^j, from W^j^^, 

oc Tr[. . . M~\t, h) '^. Q2~_^[il'ifJD2 ^7.M"^(t2, 0) ■ ■ ■], (14) 

where we have used r = it and Dt = iDr so that D = {D, —iDt) is the covari- 
ant derivative in Minkowski space. Expanding T^y in the unphysical region where 
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< 1, the expression between the 7's in Eq. (|T^) gives 
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Figure 3: Quark skeleton diagrams in the Euchdean path integral formalism consid- 
ered in the evaluation of matrix elements for the sum of local operators from the 
operator product expansion of J^(x) J;^(0). (a), (b) and (c) corresponds to the oper- 
ator product expansion from Fig. 1(a), 1(b) and 1(c) respectively. 

From this we obtain the valence and CS parton leading twist contributions to T^y 
from Fig. 1(a) 

T,.(g.^..) = E ^)\^v,v. E ^'l^^yT ^W.i) - 2 V E ^^^A;;(a/.)] + ... 

(16) 
where / indicates flavor. For the nucleon, it only involves u and d in the C I. A^IC.I.) 
is defined through the following consideration. We first note that the short-distance 
expansion in Eq. (|15D leads T^u{qv+cs) to a series of terms represented by the three- 
point functions in Fig. 3(a) which correspond to matrix elements calculated through 
the C. I. expression 

d[A]detM{A)e-^<^Tr[- ■ ■ M-\t, t2)0]M-\t2, 0) ■ ■ ■], (17) 

where the operator OJ is 



0] = ^l,l{ — r-'D,2D,3...D 



fin • 



The ratio of these three-point function in Fig. 3(a) to the appropriate two-point 
functions then define the forward matrix element which in turn defines the coefficient 
A]{C.I.) in Eq.(|16|), 

{p\^0]^\p)c.i. = A]{C.I.)2p,,p,,...p,„. (19) 

Similarly, we can perform the short-distance expansion for the CS antiparton in 
Fig. 1(b) and obtain the same expression as in Eq. (|16D except with the substitution 
q — * —q. As a result, this leads to the even n terms minus the odd n terms instead of 
the sum as in Eq. (|16D, i. e. 

T,A<lcs)= E ■■■A]{C.I.)- E ■■■A-{C.I.). (20) 

even=2 odd=3 

In other words, the short-distance expansion of T^^i, from Fig. 1(b) yields three-point 
functions with a series of insertion of the same operators ^0^"$ except with a minus 
sign for the odd n terms. This is illustrated in Fig. 3(b) where F^ denotes the n-th 
term insertion with the operator OJ and the associated kinematic factors in Eq. ([T6|) . 
Comparing F^ in Fig. 3(b) and the corresponding F„ in Fig. 3(a), the minus sign for 
the odd n terms in Eq. (^) implies that F^ = (— )"'F„. 

By the same token, the short-distance expansion for the DS parton/antiparton 
contribution in T^j, from Fig. 1(c) gives 

T^.(Wto)= E ■■■A}{D.I.) +/- E ■■■^/(^•^•)- (21) 

even=2 odd=3 

They have the same expression as T^uiQv+cs) and T^uiQcs) except AJ{D.I.) are from 
the D. I. part of the matrix element 

{p\m]<il\p)n.i. = A]{D.I.)2p,,p,,...p,„. (22) 

In this case, the leading twist expansion of the DS contribution to T^^^ in Fig. 1(c) 
now leads to two series of forward matrix elements of D. I. One is for T^^lq^s) with 
even plus odd n terms; the other is T^y{qds) with even minus odd n terms as given 
in Eq. (pi]). Both are represented in the three-point functions in Fig. 3(c). It is 
worth pointing out that T^yi^q^j^cs) in Eq. (pjf ) and T^u{(lds) in Eq. (|2l|) are the same 
as those derived from the contraction of the inner pair of the quark fields in the 
conventional operator product expansion of the time-ordered current-current product 
g(x)7^g(x)g(0)7i.g(0) [|]. On the other hand, T^y^q^s) in Eq. (jg^) and T^y{qds) in Eq. 
(pip are the same as those from the contraction of the outer pair of the quark fields in 
the current-current product. The only difference is that the path-integral formalism 
allows the separation into the CI. and the D.I. 

When the parts in Eqs. (|T^), (PP]), and ( pij ) are summed up, only the even n 
terms of the OPE are left 

= 2 E ■■■iA]iC.I.)+A]iD.I.)). (23) 

n=2,even 



This is the same as derived from the ordinary OPE. However, what is achieved with 
the path-integral formulation is the separation of C. I. from D. I. in addition to the 
separation of partons from antipartons. This has not been possible with other for- 
mulations, e.g. the light-cone definition of the distribution function. This separation 
facilitates the derivation of the different small x behavior between the CS and the DS, 
the identification of the CS parton as the source of the Gottfried sum rule violation, 
and a different evolution of (jcsix^Q"^) from that of qds{x,Q^) and qds{x,Q'^) as we 
shall see. 

Now consider the contour integral of T2 in the C I. and D. I. parts of T^,y around 
1/ = in the complex u plane while keeping Q^ fixed. / ^ ^J^!_l ' picks up the i^"~^ 
term in the series expansion of T2 in Eqs. (0), {^^, and ( PT| ) (N.B. —2q ■ p = 2MpV 
in these equations.) 



j;du2MpT2{u,Q^) _y-o„2,2M, 

/ 



ra.sM,iA..^-, . V84(iiir-M5, (24) 



for both the C I. and the D. I. The contour of the integral can be distorted to turn the 
above integral into an integral over the discontinuities of T2. Through the dispersion 
relation, this gives 



00 



2 , du2Mp2iImT2{u,Q') ^ g.2M,y„-i f' n-2 '^MpuW2{u,Q') 
Jq2 2m z/"-i ^ Q^ ^ Jo 4 ■ y ) 

Equating these two integrals and relating W2 to the parton distribution function, 
one obtains the moment sum rules. Since T2(Q^,i^) is symmetric w.r.t u — >■ —ly, we 
obtain only the sum rules for even n. Thus, 

ylf--(D./.) = M"(D./.) = £ dxx''-\qdsix,Q^) + qdsix,Q^)). (26) 

Similarly we obtain the moment sum rules for the odd n from the W3 form factor 
through the interference of the vector and axial vector currents, 

A^f'^^CI.) = M"{C.I.) = jyxx^-\q,+,s{x,Q') - qcs{x,Q')), 
AY^'^DJ.) = Kr\D.L) = f^ dxx''-\qds{x,Q'') - qdsix.Q^)). (27) 

One can define the valence parton distribution 

g^(x, Q"^) = qv+cs{x, Q^) - qcs{x, Q^). (28) 
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In this case, A^ "'^ gives the sum rules for the valence distribution. In particular, 
valence number sum rules 

Ml{C.I.)= I dxu,{x,Q^) = 2, 
Jo 

Ml{C.I.) = I dxd^{x,Q^) = 1, 
Jo 

MJ{D.I.) = 1^' dx {qds{x, Q^) - Qdsix, g')) = 0, (29) 

for the u and d quarks in the proton reflect the charge conservation of the vector 
current \I^Z7^\& and the fact that the DS carries no net charge. 

We note that the matrix elements associated with A'j^'^^'^"'{C.I.) include not just 
the valence but also the CS contribution. This is why the matrix element y4^^(C./.), 
which corresponds to the momentum fraction M^{C.I.) = {x)c.i., when calculated on 
the lattice |^, |T0[, is larger than that obtained from the experiments for the valence 
partons only, i.e. {x)c.i. > {x)v at /i ~ 2 GeV. 

3.1 Operator Rescaling and Mixing, and Parton Evolution 

The dimensionless coefficients A'^ are not constants, but depend logarithmically on Q^, 
the renormalization point of the operator product expansion. The operator rescaling 
and mixing analysis [0, |1^ for the twist-two flavor non-singlet and singlet opera- 



tors gives the renormalization group equations for the corresponding moments of the 
structure functions. In the context of the present path-integral formulation of OPE, 
the non-singlet coefficients A^ and the non-singlet moments only have contributions 
from the C I. (Fig. 3(a) or Fig. 3(b)), since their D. I. contributions cancel among 
different flavors. On the other hand, the singlet coefficients A"^ and the singlet mo- 
ments have both the C. I. and D. I. (Fig. 3(c)) contributions. Therefore, it is possible 
to go to the single flavor basis and classify the equations in terms of C I. and D. I. 
For the C I. which involves only the valence flavors (e.g. u and d for the nucleon), 
the renormalization equation is 

_J_M;(a/.) = ^<,M;(a/.). (30) 

where a^ is the anomalous dimension coefficient. For the D. I. the equation is 

d a iO'^'\ 1 + i-Y 

j^^,M]{D.L) = -^{alM]{D.I.) + ^^^^Mg), (31) 

where a"^ is the anomalous dimension coefficient for operator mixing with the gluon 
operators and Mq is the moment for the gluon distribution function. We note that 
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this mixing with gluon operators only contributes to n = even. As we can see, 
Eq. (|3lD involves the DS parton only. 

We should stress that in the literature [|^, |T6|, [l^ the non-singlet case has fre- 
quently been identified with valence. This is clearly incorrect. As we see from Eq. ( ^6]) 
that A^c^^'"'^^ [C .1 .) includes the CS partons in addition to the valence partons. De- 
tailed study of this subject on the lattice has been carried out for the matrix elements 
and form factors of the nucleon from the three-point functions as well as hadron masses 
from the two-point functions [Q. It is shown when the CS quarks are removed by 
prohibiting pair-production through the Z graphs in the C I., the hadron structure 
and masses are greatly affected. It is learned that the CS quarks are responsible for 
the meson dominance in the form factors, the deviation of Fa/ Da and Fs/Ds from 
the non-relativistic SU{6) limit, the hyper-fine splittings, and the constituent quark 
masses @]. In the context of the parton model, they are responsible for the difference 
between u{x) and d{x) in the proton. 

Following Altarelli and Parisi [|l^ , the rescaling and mixing equations in Eqs. ( ^0|) 
and (pT]) can be translated into integral-differential equations which are the evolution 
equations for the parton densities. Therefore, for the C. I. the evolution equation for 
the unpolarized valence and CS partons is 



I —P,,{-)q.+cs{y,Q'), (32) 

Jx y y 



rflnQ2 2iT Jx y ^^'y 

where 

l'dzz-'P,,{z) = ^. (33) 

For the CS antiparton density, the equation is similar 

dqcs{x,Q^) as{Q^) f^dy x ^ 

dmQ^ 2iT Jx y y 

For the DS partons in the D. I. (Fig. 1(c)), the evolution equations from Eq. (^) are 

dh^Q^ = ^^Jx Y^P^^(-)(^^s + qds){y,Q) + P,c{-)G{y,Q )], 

(35) 
d{qds-qds)ix,Q'^) a,(Q2) f^ dy x _ ^^ /q^n 
dh[Q^ = ^^Jx y^-(-)(^--^-)(^'^ )' (36) 



where 



^ PqG[Z) - — 

and G{y, Q"^) is the unpolarized gluon distribution function. 



dzz"-^P,G(^) = ^, (37) 

4 
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Finally, the gluon evolution equation is 



d In Q"^ 2tt Jx y y 



•ix y y J. , r, j-_noj-7. 



f=val.fla. f=DSfla. 



,X^ ^, 2\ 



+ PGG{-)G{y,Q')}, (3^ 

y 

It appears that, except for the gluon distribution, the evolution equations derived 
above are different from the DGLAP evolution equations a la Dokshitzer [^, Gribov 



and Lipatov |T^, and Altarelli and Parisi [|T3l due to an extra CS degree of freedom 



in Eq.(|3^ which evolves like the valence in Eq. (p2D. However, since the evolution 
equations are linear in the parton distribution functions, once one defines the total 
sea to be the sum of CS and DS, i.e. 



qs{x, Q^) = qcs{x, Q"^) + qds{x, Q"^) 
Qs{x, Q^) = qcs{x, Q"^) + qds{x, Q"^) 



(39) 



the sum of twice of Eq. (p^ and Eq. (|35D becomes 

,1 ^2 = ^ / —\Pm{-ms + qs)[y.Q) + PqG{.-)G[y,Q)\. 40 

drnQ-^ 271 Jx y y y 

This is exactly the evolution of sea partons in the DGLAP equation. Similarly, the 
difference between Eq. ( p2D and Eq. (|34D leads to the evolution of the valence partons 
in the DGLAP equation, 

dqy{x,Q'^) as{Q^) /■! dy 



f ^PU-Kiy,Q')- (41) 

Jx V V 



rflnQ2 271 Jx y '^'y 

Thus, we see DGLAP equations can be derived from the path-integral formalism. 
However, the present form with the separation of CS from the DS offers a separate 
evolution of the CS in Eq. (0) which are decoupled from the valence, the DS, 
and the gluon. This decoupling is retained when higher orders are considered ||18|| . 
This affords the possibility of maintaining the separation of the CS and DS parton 
distributions when they are fitted to experiments at certain Qq and evolved up or 
down in Q^. This is consistent with the classification into valence, CS, and DS in 
Fig. 1 and Eqs. (|^) and (^) at all Q^ from the path-integral formalism in the first 
place. Although, at leading twists, it makes no difference whether one considers the 



separate evolutions of CS and DS in Eqs. (|3^ and ( pSQ or the combined evolution in 



Eq. (^OD, it may be necessary to consider the separate evolutions of CS and DS when 



higher twists are taken into account [jT9|. 
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4 Phenomenological Consequences 

After the dynamical parton degrees of freedom are classified as valence, connected 
sea (CS) , and disconnected sea (DS) via the path-integral diagrams in Fig. 1, one 
may question if there is an advantage to separating the CS from the DS. After all, 
both the CS and DS are pair-produced sea, only with different topology in the path- 
integral diagrams. In the previous section, we have shown that the CS parton evolves 
the same way as the valence, i.e. it is decoupled from the DS and the gluons in the 
evolution equation. This is consistent with the classification of the parton degrees of 
freedom into valence, CS, and DS at all Q^. In the following, we shall show that it is 
useful to distinguish CS from DS because they have different small x behaviors. This 
is especially important in view of the fact that, in the nucleon, u and d have both 
the CS and DS parts, yet s has only the DS part. We further note that the Gottfried 
sum rule violation is attributable to the CS partons primarily. We finally asses the 
magnitude of the momentum fraction due to the CS and DS partons. 

4.1 Small X behavior 

Since the CS parton is in the connected insertion which is flavor non-singlet like the 
valence, its small x behavior reflects the leading Reggeon exchanges of p, a;, 02-.. and 
thus should be like x~^/^ pO, pTI, WA p3|. On the other hand, the DS is flavor singlet 



and can have Pomeron exchanges, its small x behavior goes like x ^ ||20|, |2l|, ^ 
Therefore, we have 

Ues/4s(a;,Qo),.-oa;"^/^ (42) 

Uds{x, Ql) ~ dds{x, Ql) ~ s(x, Ql)^z,^x~'^- (43) 



at certain Qq. As a result. 



u{x,Ql)-d{x,Ql)^^^x-^/\ (44) 



and this has been incorporated in the fitting of experimental results at some input 
scale Ql [^ ^ ^. On the other hand, it has been taken for granted that u{x, Ql) + 



(i(x, Ql) has the same behavior as s(a;, Qo) in the fitting of the experiments. In other 
words, it is assumed that 

^(m(x, Ql) + d(x, Ql)) ~ 2s(x, Ql), (45) 

in the fitting of parton distributions at Ql = 1.6 GeV^ for CTEQ4 |Q and Q, 



1 GeV^ for MRST [Q. This has made the assumption that both the u, d, like s, have 
only the DS component. We proved in Section ^ (Eqs. (|^) and (^)) that this is not true 
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and that u and d have CS in addition to DS. Therefore, the correct parametrization 
for |(u(x, Ql) + d{x, QD) should be 

-(n(x, Ql) + d{x, QD) = As{x, Ql) + -{u,s{x, Ql) + 4s(x, Q^)), (46) 

where A is a proportional constant. In the SU{3) flavor symmetric limit, A = l.The 
second term on the right hand side of Eq. (|4^) is the CS part and has the x~^^'^ small 
X behavior. As far as we know, this kind of parametrization for ^(^(x, Qq) + J(x, Qq)) 
has not been taken into account in extracting the parton distribution functions. When 
this CS degree of freedom is incorporated in the sum of u and d in addition to the 
difference, we will have a different result from the present global analyses which will 
lead to different predictions on the parton distributions at high Q^ relevant to LHC. 



4.2 Origin of Gottfried sum rule violation 

The NMC experiments of the F2 structure functions of the proton and deuteron ^ 
reveal that the Gottfried sum rule 

Sa=r,.^M^I£^^l (47) 

JO x 3 

is violated due to fact that u y^ d in the proton. This is verified in the Drell-Yan 
experiment E866/NuSea [0. It has been shown that the DS partons in Fig. 1(c) 
cannot give rise to a different u and d when m„ = m^. It is noted that in Fig. 1(c) the 
flavor indices in the quark loop are separately traced from those in the propagators 
associated with the nucleon interpolation field, only the latter reflects the valence 
nature of the proton. Hence, Fig. 1(c) does not distinguish a loop with an u quark 
from that with a d quark at the flavor symmetric limit, i.e. m„ = rud- In other 
words, Udsix, Q"^) = dds{x, Q^). On the other hand, the origin of u{x, Q^) 7^ d{x, Q^) 
can come from the CS antipartons in Fig. 1(b). Thus, the violation of the Gottfried 
sum rule originates entirely from the CS partons in the charge symmetric limit, i.e. 

X 6 6 J 

We shall see later that to 0{as) and 0(1/ Q"^), this turns out to be a sum rule. 

Due to the fact that 

dxP,,{x) = ^ = 0, (49) 



/ 

JO 



where aj^ is the anomalous dimension coefficient in Eq. (PH]), Eqs. (^TJ) and ( [3^ ) lead 

to 

d /■! 

dxq^,{x,Q^) =0, (50) 

rfxfe(x,g2) = o. (51) 



dlnQ^ Jo 
d /■! 



din Q^ Jo 
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We see that the CS antiparton number, hke the valence number, is conserved, i.e. 
independent of Q^. This is in contrast with the DS partons whose number is not 
conserved due to the pair-creation from the gluon. However, the conservation of 
the CS antiparton number is only good in the leading logarithmic approximation, i.e. 
good to 0{as) and 0{1/Q'^). Whereas, the conservation of the vector current protects 
the charges, thus the valence quark number, against any Q^ correction [|ll|, [1^, [13|. 

We note that the sum in Eq. ( P5| ) is in terms of the CS antipartons numbers. 
Thus to leading logarithmic approximation, it is a sum rule 

'^^ \. = 3 + 3 K. - ridj. (52) 

where riu^Jn^^^ is the Ucs/dcs number. 

4.3 Magnitude of the connected sea partons 

We don't know precisely how large the magnitude of the CS partons and antipartons 
are in comparison with those of the valence and the DS unless one fits the DIS and 
Drell-Yan experiments with an explicit separation of the CS and the DS. But there 
are hints which are helpful in this respect. From the NMC experiments on the F2 
structure function of the proton and deuteron 0, the sum in Eq. (|52D is measured 
to be 0.235 ± 0.026 at Q^ = 4 GeV^, significantly smaller than the Gottfried sum rule 
prediction of 0.333. This implies that 

nu.. - n^, = / dx[ucs{x, Q^) - d,s{x, Q^)] = -0.147 ± 0.039, (53) 

at Q^ = 4 GeV^ which is not negligible compared with the valence numbers of u and 
d. 

Since it (similarly d) has contributions from the CS and DS, i. e. u{x) = Ucs{x) + Uds{x) 
and s is from the DS only, one expects that 

{x)u = {x)u,, + {x)u^^ > {x)-s. (54) 

Indeed, in the CCFR dimuon data, s is about 50% of {u + d)/2 at Q^ — 4 GeV^ in the 
range of x between 0.01 and 0.20. Since u — d, which only has the CS contribution, 
also peaks in this x range, we think the observed difference between {u + d)/2 and s is 
mainly duw to the CS part of {u + d)/2. As a result, one expects that the momentum 
fraction of the CS part of (n + d)/2 is comparable to that of the DS part, i.e. 

at Q^ = 4GeV^. Any experiment which measures {u + d)/2 and s at very small x, 
e.g. 10"'^ — 10~^ will be very useful in verifying the form of the distribution function 
prescribed in Eq. (^61). 
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5 Conclusion 

In conclusion, we have formulated the hadronic tensor W^,y of the deep inelastic 
scattering starting from the Euclidean path-integral formalism. We found that it 
can be divided into three gauge-invariant and topologically distinct parts which we 
classify as the valence-connected sea partons, the connected sea antipartons and the 
disconnected-sea partons and antipartons. This admits a separation of the C I. from 
the D. I. and the partons from the antipartons. Since the CS is in the C I. and the DS 
in the D. I., they have different small x behaviors. We show that the operator product 
expansion is simply a short distance Taylor expansion of functions in the path-integral. 
From operator rescaling and mixing, we derive the evolution equations which show 
that the CS partons evolve like the valence and their numbers are conserved in the 
leading log approximation. We stress that in the nucleon u and d partons have both 
the connected and disconnected sea contributions; whereas, s partons has only the 
disconnected sea contribution. A global analysis of the experimental data is needed 
to take this difference into account to fit the parton distribution functions. 
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L. Mankiewicz, R. McKeown, CS. Lam, G. Martinelli, J.C. Peng, J. Qiu, W.K. Tung 
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